Prom band insulator to Mott insulator in one dimension 
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We derive the phase diagram for the one-dimensional 
model of a ferroelectric perovskite recently introduced by 
Egami, Ishihara and Tachiki [Science, 261, 1307 (1993)]. We 
show that the interplay between covalency, ionicity and strong 
correlations results in a spontaneously dimerized phase which 
separates the weak-coupling band insulator from the strong- 
coupling Mott insulator. The transition from the band insu- 
lator to the dimerized phase is identified as an Ising critical 
point. The charge gap vanishes at this single point with the 
optical conductivity diverging as a-{uj) ~ o;^'''^^. The spin ex- 
citations are gapless above the second transition to the Mott 
insulator phase. 
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Much effort have been devoted to understanding the 
response of strongly correlated electron systems to lat- 
tice distortions. Generally, one expects that the tendency 
to a charge density wave instability modulated by a lat- 
tice distortion is increased by (repulsive) interactions. A 
surprisingly pronounced enhancement has recently been 
observed in numerical studies of a prototype one dimen- 
sional (ID) model for ferroelectric perovskites, proposed 
by Egami, Ishihara and Tachiki [|]. 

The Hamiltonian of their model is 

^ = E + H.c.) -I- A(-l)'n,,^ 

+ U^n^in.,^l, (1) 

i 

at half-filling. The odd and even sites represent oxygen 
atoms (O) and a generic cation (C), respectively, with the 
energy difference Ec — Eq — 2 A. The hopping ampli- 
tude t determines the amount of covalency and U stands 
for the on-site Coulomb repulsion. In spite of its ap- 
parent simplicity, model (^ reveals nontrivial physics. 
At ?7 = it describes a band insulator (BI) of a mixed 
ionic-covalent character, with a spectral gap for all ex- 
citations. At A = it is a Mott insulator (MI) with a 
finite charge gap, rric, but with gapless spin excitations. 
In fact, in the strong-coupling limit U ^ (^, A), (^ can 
be mapped onto a Heisenberg model by projecting out 
all doubly occupied sites. At the leading order, only a 
nearest neighbour exchange J = 4t^L//(C/'^ — 4A^) is gen- 
erated. The resulting Heisenberg chain is well known to 
possess a gapless spin spectrum. The issue of interest is 



the nature of the crossover from the BI regime to the MI 
regime which, on general grounds, is expected to occur 
in the region where the single-particle BI gap A becomes 
comparable with the MI charge gap rUc. 

Thus, increasing U at fixed A, one expects a transi- 
tion to a spin-gapless phase. In addition to this spin 
transition, the charge degrees of freedom should be also 
affected in the course of the crossover. This can easily be 
understood in the case t <C {U, A) where, upon increas- 
ing U , the system approaches a mixed-valence regime 
in which the two charge configurations, 0~^C"'"^ with 
energy U — 2A and 0~"'^C+^ with zero energy, become 
degenerate. Therefore, those excitations responsible for 
the charge redistribution among these two configurations 
should soften around some value of U . The finite-size 
simulations suggest that the mixed-valence regime 
is accompanied by a strongly enhanced response to the 
coupling with zone-centre optical phonons ^, 
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H.C, 
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where A is the electron-phonon coupling constant and 
Vi is the dimerization (or staggered bond-density) op- 
erator. In fact, the so-called average dynamical charge 
= lim^^o{P{^) - P(0))/^, where P{^) is the polar- 
ization, diverges at a particular point, even in finite-size 
simulations |4|^. This is due to the accidental degener- 
acy between two singlet ground states with opposite par- 
ity which occurs in any finite-size system with the num- 
ber of sites multiple of 4 [^,||. More importantly, Resta 
and Sorella have recently observed that the localiza- 
tion length of the ground-state wave function diverges at 
some value of U, which suggests vanishing of the charge 
gap; hence the charge transition. 

However, the overall picture which emerges from nu- 
merical simulations and the existing mean-field calcula- 
tions 1^ , clearly underestimating the crucial role of fiuc- 
tuations in ID, is still far from being satisfactory. It 
remains unclear whether the charge instability occurs at 
the same point as the spin transition. Even the nature of 
the transition (transitions?) is not known. Is it a first- 
order transition or is it a second-order one? In the latter 
case, what is the universality class it belongs to? To 
the best of our knowledge, there have been no consistent 
analytic attempts to resolve these intriguing questions. 
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In this Letter we present theoretical arguments show- 
ing that there are two continuous transitions: the spin 
transition of the Kosterlitz-Thouless (KT) type at [/ = 
Uc2, and the charge transition at J7 = Uci < Uc2, the lat- 
ter identified as an Ising critical point where the charge 
gap vanishes. In the intermediate region, Uci < U < Uc2, 
the site-parity is spontaneously broken, and the system is 
characterised by a doubly degenerate, dimerized ground 
state, with {Vi) being the order parameter. 

We start our discussion with addressing the mecha- 
nism of the spin gap generation in the MI phase. For the 
single-chain Hamiltonian (|l|), which is SU(2) and site- 
parity invariant, the only possibility compatible with the 
symmetries is a spontaneous spin-dimerization, like the 
one occurring at J2/J1 — 0.24 in the frustrated Heisen- 
berg chain with nearest-neighbour (Ji) and next-nearest- 
neighbour (J2) exchanges. The same scenario is realized 
in our case below a critical value Uc2 which depends on 
A and t. This result agrees with the exact diagonaliza- 
tion of a 4-site chain where the next-to-nearest neighbour 
exchange starts to play a role. At large C/, the ground 
state is an odd-parity singlet, and the lowest excitation is 
a spin-triplet, followed at higher energy by an even-parity 
singlet. As U decreases, first the even-parity singlet and 
the triplet cross, and finally the two opposite-parity sin- 
glets cross. For ([/, A) ^ t these two level crossings oc- 
cur at [/ ~ 2 A 1.86 1 and at f7 ~ 2 A + t, respectively 
(see also Ref. Q). The same is observed in the J1-J2 
Heisenberg model when J2 is increased 0. The crossing 
of the triplet and singlet excited states in finite-size sys- 
tems becomes the spontaneous spin-dimerization transi- 
tion in the thermodynamic limit. The crossing of the 
two opposite-parity singlets, taking place at J2 = 0.5 Ji, 
corresponds to the Majumdar-Gosh solvable limit 0. 

The appearance of a spontaneously dimerized insulat- 
ing (SDI) phase can be also inferred starting from the BI 
phase. In this limit we need an improved version of the 
mean- field calculation presented in Ref. We intro- 
duce a two component Fermi field ^fj^^ = (c^^ , c^^^^), 
fc G [— 7r/2, TT, 2], and the operators 

TaM-^Y.'^L^-'i'u+.a^ (a -0,1, 2, 3), (3) 
fc 

where Tq are the Pauli matrices including the unit matrix 
tq. The staggered charge/spin densities correspond to 
fif (0)±fi|(0), while the dimerization is 7^21 (0)-|-f2x(0) = 
[seeEq.(^]. A generic interaction compatible with 
the symmetry of (|l|) can be written as 

Htnt ^ jj- [(51 + 93)^i<T{q)n<7'{-q) 

+ {91 - 93)T2a{q)T2<T'{~q) + 2g2To^{q)TQ^'{-q)]. (4) 

The coupling constants are written using the (7-ology jar- 
gon j|]: gi,g2 and 53 label the backward, forward and 



Umklapp scattering amplitudes, respectively. In the spe- 
cial case of (0), 31 = 52 = ff3 = U, and the f2-f2 term 
vanishes. Therefore, in any mean-field treatment of the 
original Hamiltonian, only variational states with finite 
average values of ricr(O), i.e. the staggered charge and 
spin densities, can be explored. However, the equality 
among the 5's is already lost at the first order of per- 
turbation theory. More refined weak-coupling schemes, 
such as renormalization group or the bosonization 
technique pO| , show that renormalization occurring at 
energies well below the bandwidth cutoff but still larger 
than ■max{A,mc) enhances and 32 but suppresses gi. 
This implies generation of an effective coupling {g^ — gi) 
which may lead to the appearance of a finite average 
value of f2(0), in addition to the average (fi(0)) induced 
by the A-term. By including these renormalization ef- 
fects within the leading logarithmic approximation 
we indeed find a sequence of two transitions. First, at 
U ~ Uci, a finite average (t2(0)) continuously appears. 
At a larger U = Uc2 a triplet exciton gets soft signalling 
the transition to the MI phase. Although these two tran- 
sitions are physically distinct within this scheme, to es- 
timate the difference between Uci and Uc2 one has to go 
beyond the leading logarithmic approximation, which is 
not feasible. For this reason, we can only safely state that 
Uci < Uc2, with both critical points being of the order 
27rt/ln(</A). 

In what follows, we develop the low- energy effective 
field theory for the lattice model (|l]). Considering the 
weak-coupling case, [U, A) ^ t, we linearise the spec- 
trum and pass to the continuum limit by substituting 
Oq ^^^c„o- i^i'Raix) + {-iY^La{x), X = noo, where 
flo is the lattice spacing, and i^R.Lix) are the right 
and left components of the Fermi field. These fields 
can be bosonized in a standard way |l^: ipR,L;(T = 
(27rao)"^/^e=^'^'^«-^'<', where are the right(left)- 

moving Bose fields. We define $0- — (j)Ra + <pLa and intro- 
duce linear combinations, <I>c = -I- ^^1)^/2 and $s = 
(<i>l — $x)\/2, to describe the charge and spin degrees of 
freedom, respectively. Then the Hamiltonian density of 
the bosonized model is given by: Ti-eff — Ti-c + 'Hs + 'Hcs- 
Here the charge and spin sectors are described by 

= I [n^ + [d.'^cf] 

cos \/87r'I>c + -^^dx<j)cRdx<l>cL (5) 

+ cos a/Stt^s —dxCjisRdxtpsL (6) 

where H^.s are the momenta conjugate to ^c.s, and Vc^s 
are the velocities of the charge and spin excitations. Ac- 
tually, Tis in is an Abelian bosonized version of the 
SU(2)i-symmetric Wess-Zumino-Novikov-Witten model 
with a marginally irrelevant (51 > 0) current-current per- 
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turbation, —giJn -Jl, given by the last two terms in (||). 
The two sectors of the theory are coupled by the A-term: 



■He 



- (2 A /TTflo ) sin V 27r$c cos v 27r$s . 



(7) 



In fact, the effective Hamiltonian li-eg represents the low- 
energy limit for a wide class of lattice models including 
(|l|). Such lattice effects as nearest-neighbour repulsion, 
or a difference between the Coulomb repulsion on the O- 
atoms and the C-atoms, will only renormalize the param- 
eters of 'H.eff^ in particular the charge stiffness constant 
(Kcj and the velocities Wc,s- 

Starting from the MI phase, one can firmly justify the 
opening of a spin gap and spontaneous dimerization upon 
decreasing U . Assuming [/ ^ A, in which case the spin 
sector is gapless while the charge one is gapped, one can 
integrate out the charge degrees of freedom to obtain at 
the second order in A a nontrivial renormalization of g\ : 
gi = gi — (A/toc)^ (vc/ao) . Until gi > 0, the spin spec- 
trum remains gapless. However, when U is decreased, gi 
eventually becomes negative, and the system undergoes a 
continuous (KT) transition to the dimerized phase with a 
nonzero spin gap. Using the exact result for the small-?/ 
Hubbard model, rric ~ \/Tlie~'^^*^^ , one concludes that 



27rt 



c2 



ln(t/A) 



1 + 



lnln(</A)\ 
In(VA) ) 



(8) 



in agreement with the previous estimate. 

Bosonization also offers a simple understanding of the 
region U < Uc2- In what follows, we shall focus on the 
nature of the transition at [/ = Uci- To this end, the 
Hamiltonian H-eff, given by will be regarded as 

a phenomenological Landau-Ginzburg energy functional, 
in the sense that all the couplings are effective ones ob- 
tained by integrating out high-energy degrees of freedom. 

We start with looking at the saddle points of the ef- 
fective potential in Ti-eff which obey the equations (here 
ifics = \/27r$c,s) 



cos ipc 



i sin ipc — 27rA cos ips 



= 



sin ifs 27r A sin tpc — ^gi cos = 



(9) 
(10) 



supplemented by stability conditions. At 53 < 7rA/2 one 
finds two sets of minima (defined modulo 27r) located at 
ifs — 0,ifc — 7''/2 , and ips — "^c — — vr/2. These sets 
characterise the BI phase. Indeed, the vacuum-vacuum 
transitions, A(y5s(c) = iTr, describe stable topological ex- 
citations carrying the charge Q = l^ipc/'K = ±1 and spin 
S'^ = ^ipsjl-K = ±1/2 and therefore coinciding with 
"massive" single-fermion excitations of the BI. The sit- 
uation changes at g^ > 7rA/2 where each minimum in 
the charge sector splits into two degenerate minima and 
thus transforms to a local double-well potential. The new 
minima are given by ips = 0, (pc = 4>o, tt— ^o, and ips — tt, 
tpc = —00, — + 4>0i where 0o = arcsin(7r A/253). They 



describe the SDI phase. In this region, the dimerization 
operator 'D(x) ^ cos tpc(x) cos i^s(x) acquires a finite ex- 
pectation value. Notice that the location of the minima 
in the spin sector, and hence the spin quantum num- 
bers of the topological excitations, are the same as in the 
BI phase (53 < 7rA/2). Therefore the spin part of the 
spectrum formed in the SDI phase smoothly transforms 
into that of the BI phase. The transition at ^3 = 7rA/2 
mainly involves the charge degrees of freedom which un- 
dergo dramatic changes. In particular, the charge quan- 
tum numbers become fractional, depending on 4)q. The 
2'2-degeneracy of the SDI state implies the existence of 
topological kinks carrying the spin S = 1/2 and charge 
Q — ±2(/3o/2. These massive excitations, described by 
interset vacuum-vacuum transition, interpolate between 
the neutral spinous of the MI phase and single-fermion 
excitations of the BI phase. However, the double-well lo- 
cal structure of the effective potential in the charge sector 
gives rise to singlet kink excitations as well, whose mass 
and charge Q = 1 — 2(/)o/7r vanish at 53 = 7rA/2. On the 
BI side of the phase diagram, these singlet excitations ap- 
pear as excitonic bound states inside the single-particle 
gap. Precisely at 53 = 7rA/2, the effective potential be- 
comes ~ (fi'^, which is well known to describe the Ising 
universality class. 

We have checked that the SDI solution remains stable 
in the region 



(11) 



This condition can be satisfied if gi < g^, which is in- 
deed true for our phenomenological Hamiltonian. How- 
ever, separation between the transition points Ud and 
Uc2 becomes even stronger when the electron-phonon in- 
teraction, always present in realistic situations, is taken 
into account. Let us assume that the zone-centre opti- 
cal phonons are dispersionless, with a characteristic fre- 
quency luq being the largest energy scale in the problem. 
In this limit, the phonons mediate an instantaneous inter- 
action between the electrons which adds to the Hubbard 
repulsion and leads to a renormalization of the "bare" 
couphng constants: 33 ^ ff3 + 5o, 5i ^ ffi - 5o, where 
go ~ /ujQ. The transition points Uci and Uc2 will start 
moving apart upon increasing go, and if the electron- 
phonon coupling is large enough (50 > 5i)i a finite spin 
gap will be generated even at A = 0. Qualitatively 
the same situation occurs upon relaxing the condition 
of extremely quantum phonons (large ujo) at a fixed A. 
This observation explains that the two transitions can be, 
in principle, completely disentangled and, on the other 
hand, it provides an extra support to the fact that the 
charge transition separates the BI and SDI phases. 

The Ising scenario at U = Ud can be more rigorously 
proven by studying the charge degrees of freedom with 
the spin bosonic field locked at = or n/2. Such 
approach will certainly be valid in a small vicinity of 
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Uci because, as we already mentioned, the structure of 
the gapped spin sector is quaUtatively the same both in 
the BI and SDI phases. The charge sector can then be 
described by a double sine-Gordon Haniiltonian 



Vc_ 

2 



2A 



cos \/87rXc$c(a;) sm \/ 271 K c{x) , (12) 

Trao 



where, we stress once more, all the parameters are purely 
phenomenological, and part of the interaction has been 
absorbed in the rescaling parameter Kc < 1. A quan- 
tum phase transition of an Ising type in the double sine- 
Gordon model was recently discovered by Delfino and 
Mussardo . In the whole parameter region of interest 
to us, i.e. 5 > and ifc < 1, they argued that the ex- 
istence of an Ising critical point is a universal property 
of this model. This allows us to investigate this model 
in the vicinity of the point Kc — 1/2, where the prob- 
lem greatly simplifies. At this point, the Hamiltonian 
( |l^ ) can be mapped onto two coupled 2D Ising models, 
or equivalently, two coupled quantum Ising chains in a 
transverse magnetic field (a deformation of the Ashkin- 
Teller model). One of the two chains always remains 
off-critical, while the other one can pass through a crit- 
ical point by fine tuning the parameters of the model. 
The details will be presented elsewhere Here we 

just quote our main findings, which can be easily un- 
derstood once the Ising scenario is assumed. The most 
relevant perturbation at an Ising transition is the mag- 
netic field h, which, in the ordered phase, selects one of 
the two equivalent vacua. By analogy, the role of h is 
played in our case by the explicit dimerization ^, Eq.(^. 
The dimerization operator T) (which also plays the role 
of the charge polarizaton operator) is therefore propor- 
tional to the order parameter a of the quantum critical 
Ising model. The Ising mapping also allows us to identify 
the electron current operator, j ^ dtcr, and the average 
electron charge operator p ^ —dx<y- These findings are 
perfectly consistent with the general relations p ~ —dxT^ 
and dtp -\- dxj — (continuity) valid for insulators (i.e., 
in the absence of free carriers). ^From known results on 
the correlation functions at the Ising transition we con- 
clude that the susceptibility to the explicit dimerization, 
namely the real part of the polarizability x{^)i diverges 
at low frequencies as uj~'^^^, and the optical conductivity, 
cr(w) ~ uj~^/^ (no Drude peak), is 'semi- metallic'. 

In conclusion, we have identified the origin of the 
anomalous response to an electron-phonon coupling in 
the ID ferroelectric perovskite model (|l|) , which turns out 
to be far stronger than that of free electrons (A = [/ = 0), 
of a band insulator {U = 0), or of a pure Hubbard model 
(A = 0). Our findings are summarised in Fig.l, where 
the phase diagram is shown at fixed A and t as a func- 
tion of U. Between the BI {U < Ud) and MI (U > Uc2) 
phases there exists a spontaneously dimerized phase even 



in the absence of an explicit electron-phonon interaction 
[|l2|. The transition from the MI to the spontaneously 
dimerized insulator is of the KT type, accompanied by 
opening of the spin gap. The transition from the dimer- 
ized phase to the BI has been shown to be of the Ising 
type, with the charge gap closing at this single point [ p^ . 
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and t. Ac{s) stands for the charge(spin) gap, 7?i is the dimer- 
ization operator, whose average is finite in the spontaneously 
dimerized insulator but vanishes both in the band and in the 
Mott insulators. 
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